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EXPERIMENTAL  STUIiy  OF  THE  TURBULENT  BOUNDARY  LAYER 
IN  GAS  MOTION  IN  AXISYMMETRICAL  DIFFUSERS 
WITH  COOLED  WALLS 

P.  N.  Romanenko  and  A.  I.  Leont'ev 
1.  Statement  of  Problem 

The  study  of  turbulent  boundary-layer  characteristics  under 
conditions  of  significant  longitudinal  pressure  and  heat-transfer 
gradients  Is  of  great  practical  value  for  the  development  of  engi¬ 
neering  methods  to  calculate  friction  and  heat  transfer  In  modern 
Industrial  Installations  and  to  elucidate  turbulent  mixing  under 
these  conditions. 

The  published  works  on  the  Investigation  of  the  turbulent 
boundary  layer  In  the  presence  of  significant  pressure  gradients 
[1,  2,  5^  6,  12]  are  mainly  concerned  with  the  dynamic  boundary  layer. 
The  well-known  seml-emplrlcal  methods  of  calculating  an  Isothermal 
turbulent  boundary  Isiyer  In  the  region  of  significant  pressure  gradients 
permit  sufficiently  reliable  calculation  of  the  dynamic  characteristics 
of  the  turbulent  boundary  layer  and  enable  us  to  determine  the  location 
of  separation  points  of  the  turbulent  boundary  layer.  However,  the 
expansion  of  these  methods  to  the  calculation  of  a  turbulent  boundary 
layer  under  conditions  of  Intense  turbulence  presently  remains  an 
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open  question.  Methods  of  computing  the  thermal  boundary  layer  which 
are  sufficiently  expanded  In  engineering  practice  [I3]j  based  on  the 
use  of  the  analogy  between  the  processes  of  heat  and  momentum  transfer 
and  also  on  the  self- similarity  of  the  laws  of  friction  and  heat  trans¬ 
fer  In  the  region  of  significant  positive  gradients  of  pressure,  are 
still  Inapplicable.  Studies  in  which  attempts  were  made  to  calculate 
theoretically,  the  influence  of  pressure  gradients  on  the  dynamic  and 
thermal  characteristics  of  boundary  layers  are  well  known  [l4]. 

However,  the  large  number  of  poorly  founded  assumptions  and  Inexcusably 
ciaubersome  formulas  significantly  reduce  the  practical  value  of  the 
methods  of  computation  proposed  In  these  works. 

Since  In  the  near  future  we  cannot  count  on  serious  progress 
in  the  application  of  statistical  theories  to  the  study  of  anisotropic 
turbulence,  experimentation  must  be  acknowledged  as  the  most  effective 
method  of  .study. 

This  work  presents  the  results  of  studies  of  the  turbulent  boundary 
layer  In  heated  air  flow  In  diffusers  of  circular  cross-section  with 
cone  angles  of  8 '*4*  and  12®.  The  axlsymmetrlcal  diffuser  makes  it 
possible  to  eliminate  the  effect  on  calculations  of  the  local  resis¬ 
tance  coefficient  and  of  other  flow  characteristics  of  a  three-dimen¬ 
sional  gas  flow  which,  as  was  demonstrated  by  P.  Clauser  [3]j  takes 
place  in  flat  diffusers.  The  use  of  experimental  data  obtained  for  a 
flat  diffuser  hampers  the  calculation  of  the  resistance  coefficient 
with  respect  to  the  Integral  relation  of  momenta.  In  addition,  the 
study  of  the  axlsymmetrlcal  flow  gradient  Is  of  great  practical  Interest. 

2.  Experimental  Installation 

A  general  schematic  of  the  experimental  Installation  Is  shown  In 
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Pig.  1.  It  consists  of  an  open  wind  tunnel  of  periodic  action  with 
detachable  nozzles  and  a  test  section.  In  addition^  the  apparatus 
Includes:  a  2R-V220  air  compressor  which  produces  l80  mVhr  (NTP) 
(final  air  pressure  of  220  atm  (abs),  an  80-tank  receiver,  a  300  Kw 
electric  heater,  an  arrangement  for  the  enclosure  and  feeding  of  ejected 
air  to  the  chamber,  a  system  for  water  cooling  the  test  section,  and 
measuring  devices. 


Pig.  1.  Schematic  of  experimental  Installation. 

Plgure  2  shows  a  longitudinal  cross  section  of  the  experimental 
device,  an  axlsjnnmetrlcal  diffuser  with  a  4 ''2'  cone  half-angle.  The 
profiles  In  which  static  and  total  pressures  as  well  as  temperatures 
of  the  flow  were  measured  are  .marked.  A  structural  cross  section  of  the 
sleeve  through  which  the  Pitot  microtube  and  thermocouples  are  Intro¬ 
duced  Into  the  channel  Is  shown. 

The  diffuser  Is  composed  of  five  Insulated  sections  with  double 
walls  which  form  a  space  through  which  cooling  water  Is  passed.  The 
water  flow  rate  Is  determined  by  the  gravimetric  method.  The  Internal 
wall  Is  made  of  2  mm- thick  copper  and  the  exterior  wall  of  2.5  nun 
steel. 
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Pig.  2.  Longitudinal  cross  section  of  the 
test  section. 

The  compressed  air  Is  forced  from  the  compressor  Into  the  receiver 
whose  tanks  are  Joined  by  prefabricated  ramps  Into  six  groups.  Each 
group  of  tanks  Is  connected  to  a  high-pressure  delivery  collector 
which  la  connected  to  a  control  valve.  The  control  valve  Is  connected 
on  the  low-pressure  side  to  the  duct  through  which  the  air  passes 
Into  the  heater.  The  required  air  pressure  In  the  heater  and  in  front 
of  the  nozzle  Is  maintained  by  the  control  valve. 

A  measuring  dlaphram  was  set  up  In  front  of  the  heater  In  order 
to  determine  the  flow  rate  of  air. 

3.  Measured  Quantities  and  Measuring  Devices 

The  following  quantities  were  measured  during  the  experiment: 
air  pressure  In  front  of  the  electrical  heater  and  In  front  of  the 
working  nozzle;  air  temperature  at  the  heater  outlet;  static  and 
total  pressure,  temperature  in  the  profiles  of  each  section  located 
at  a  distance  of  30  mm  from  the  plane  of  the  Joint  of  the  sections; 
static  pressure  in  two  Intermediate  profiles  of  the  first  section  In 
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the  direction  of  air  motion  in  which  the  pressure  gradient  is  greatest; 
temperatures  of  the  inner  walls  of  the  diffuser  in  the  same  profiles 
in  which  velocity,  temperature,  and  static  pressure  of  air  flow  are 
measured;  temperature  difference  and  water  flow  rate  in  all  sections 
of  the  diffuser;  pressure  difference  in  the  measuring  diaphragm. 

The  air  pressure  In  front  of  the  heater  and  working  nozzle  was 
measured  with  a  Burdon  manometer  having  a  measurement  limit  of  50  atm, 
division  value  of  0.25  atm  and  accuracy  grade  of  1.5*  The  pressure 
difference  In  the  measuring  diaphragm  was  measured  with  a  DT-150 
differential  manometer.  Total  pressure  in  the  boundary  layer  In  the 
unperturbed  flow  In  the  sections  In  front  of  the  cooled  sections  was 
measured  with  total-pressure  pneuniometric  microtubes.  One  end  of 
the  Pitot  tube  faces  toward  the  Incident  flow  and  the  other  Is  connected, 
to  a  differential  manometer  which  is  made  In  the  form  of  a  U-shaped 
tube  filled  with  ethanol  or  water. 

Static  pressure  was  measured  with  a  static-pressure  sampling 
tube  Introduced  to  the  Inside  of  the  passage  through  the  wall. 

This  tube  Is  located  In  the  same  plane  of  the  channel  cross  section 
In  which  the  open  end  of  the  Pitot  tube  Is  located.  The  end  of  the 
static-pressure  sampling  tube  Is  connected  to  a  differential  manometer. 
The  second  tube  of  this  manometer  Is  open  to  the  atmosphere. 

In  order  to  determine  the  dynamic  pressure,  the  Pitot  tube  and 
the  static-pressure  tube  are  connected  to  a  differential  manometer 
which  records  the  difference  of  levels  of  the  working  fluids  on  Its 
scale. 

The  Pitot  tubes  and  thermocouples  are  set  up  by  means  of 
traverse  gears  which  made  It  possible  to  measure  the  dynamic  pressures 
and  temperatures  In  the  control  profiles  of  the  sections  of  the  test 
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zone  at  0.05  mm  Intervals.  It  is  possible  to  measure  the  quantities 
mentioned  at  the  corresponding  points  of  all  the  five  sections  at  the 
same  time. 

The  flow  temperature  was  measured  by  chromel-alumel  thermocouples 
Installed  at  a  given  distance  from  the  wall  by  the  traverse  gear. 

Air  temperature  was  recorded  automatically  on  the  tape  of  an  electronic 
tape  balancing  potentiometer.  The  potentiometer  was  connected  to  10 
thermocouples,  5  of  which  measured  the  air-flow  temperature  and  5 
the  heat  of  the  cooling  water  In  the  sectional  housings. 

The  temperatures  of  the  Inner  walls  of  the  sections  were 
measured  with  chromel-alumel  thermocouples  placed  In  machined  grooves 
and  carefully  embedded  flush  to  the  air  flow.  Thermocouple  readings 
were  recorded  with  a  potentiometer. 

The  temperature  at  the  heater  outlet  was  monitored  by  a  chromel- 
alumel  thermocouple  connected  to  an  electronic  potentiometer  type  EPD. 

ETA-5A  electromanometers  were  also  used  to  measure  velocity, 
direction,  and  temperature  of  air  flow. 

4.  Method  of  Treating  Experimental  Data. 


In  order  to  determine  the  variation  of  local  friction  values  on 

the  wall  T  and  unit  heat  flow  at  the  wall  q,,  over  the  length  of  the 
w  w 

diffuser,  we  used  the  Integral  relationships  between  momentum  and 
energy  which  for  axlsymmetrlcal  gas  flow  Is  written  In  the  following 
form. 


For  momentum  transfer 

^  4-  , 

\  r  dx"^ 


Pi 


dx-} 


isu 

2  • 

Plttl 


(1) 
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(2) 


For  energy  transfer 

dip  r\  du,  I  dr  1  ,  d(T,,-rj, 

.rf.v"''L«i  dx'^  r  dx'^Toi  —  f,  dx 

_ _ . 

^pPi^hiTn  —  f»)  ’ 

< 

Where  r_Pfi_ /•  _^coisp\  thickness  of  momentum  loss; 

Jpi«iV  ‘hl(  .  r  j  ^ 

“  » 

pu  \  I ^  ycospN  ^  displacement  thickness; 

Pi  «i  j  \  r  j  ^ 

u  =  velocity  In  the  direction  of  the  x-axls; 

Ui=  velocity  at  outer  edge  of  the  boundary  layer; 

X  *  distance,  parallel  to  the  wall; 

r  =  radius  ( distance  from  the  surface  of  axlsymmetrlcal  channel 
to  Its  axis)  ; 

p  1=  static  density  at  the  outer  edge  of  boundary  layer; 

p  =  static  density  at  a  given  point  of  boundary  layer; 

T  =  tangential  stress  at  wall; 
w 

q  =  specific  heat  flow  at  wall; 
w 

c  «mass  thermal  capacity  at  constant  pressure; 

P* 


■  + 


Pi  dx\ 


To  =adlabatlc  stagnation  temperature  In  unperturbed  flow  ; 

T^=  wall  temperature; 
a 

('  pu  /  fo -- T’o  \ >'cosp\  , 

Pi~^  (foi  —  f  /  \ - r — thickness  of  energy  loss  or  thlck- 

b 

ness  of  enthalpy  loss; 

To=  stagnation  temperature  at  a  given  point  In  the  boundary  layer; 
To  i*equlllbrlum  temperature; 

P  =  angle  between  the  normal  to  the  wall  and  the  axis  of  the 
diffuser  In  the  plane  of  the  control  profile; 

6  ■  thickness  of  the  dynamic  boundary  layer; 

A  ■  thickness  of  the  thermal  boundary  layer. 

All  of  the  quantities  which  enter  into  Eqs.  ( l)  and  (2)  with  the 
exception  of  t  and  q„  were  determined  from  experimental  data  for 

W  Yi 
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constant  values  =  O.0705  and  O.IO5. 

Velocity  profiles  and  actual  flow  temperature  profiles  were  cal¬ 
culated  according  to  the  measured  temperature  and  dynamic  pressure 
profiles. 

Velocity  was  calculated  by  the  formula 


where j  for  air. 


u  =  d*k. 


a*  =  18,3  /To; 


True  air  temperature  was  determined  by  the  equation 


where  =  equilibrium  temperature  of  the  thermocouple 
r  =  temperature  recovery  factor; 

I  =  mechanical  equivalent  of  heat. 

A.  A.  Gukhman,  N.  V.  Ilyukhin  et  al.  [4]  studied  the  dependence 
of  the  recovery  factor  r  on  Reynolds  and  Mach  numbers  M  for  a 
longitudinally  stresimllned  thermocouple  d  *  0.5  mm  placed  In  a  tube 
with  a  diameter  of  25  mm  in  the  presence  of  a  boundary  layer.  It  was 
established  that  over  the  Interval  5.8*  10®<  Re*  <  l4. 4. 10®  and  0.24< 
<M<0.98  the  quantity  r  Is  a  constant  and  equal  to  O.9O. 

In  our  tests,  longitudinally  streamlined  thermocouples  with  a 
wire  diameter  of  0.5  mm  and  junction  diameter  of  I.5  nrni  were  used. 


The  values  of  Re* 


(where  d  Is  the  wire  diameter)  varied  from 


900  to  10,000.  Therefore,  we  assumed  r  =  0.90. 

Prom  the  measured  temperature  and  velocity  profiles  we  determined 
the  values  of  conditional  boundary  layer  thicknesses  6,  6*,  and 
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9  of  other  quantities  which  enter  Into  the  Integral  relationships  of 
momentum  and  energy  (ui,  Pi,  (Tci-Tj. 

Then  the  tangential  stresses  at  the  wall  and  the  specific  heat 
flows  were  calculated  by  Eqs.  ( l)  and  (2)  . 

« 

The  reliability  of  determining  the  local  values  of  and 
by  the  described  method  depends  substantially  on  the  accuracy  of  the 
graphic  differentiation  of  the  original  parameters.  Therefore,  local 
values  of  (or  c^  and  q^  were  calculated  by  two  Independent  methods. 
The  method  of  determining  c^  is  based  on  the  well-known  fact  that  the 
Interior  of  the  turbulent  core  of  the  boundary  layer  contains  a 
xinlversal  logarithmic  velocity  profile  which.  In  the  presence  of 
Isothermic  flow,  is  described  by  the  equation 


^  =  5,751q^*  +  5,5.  „  (7) 

where  v^^  ■  'the  dynamic  velocity  at  the  wall;  and 

V  ■  kinematic  viscosity. 

As  Is  known,  law  (7)  was  obtained  by  Integration  of  Prandtl’s 
equation  for  tangential  stress 

■  (8)  [ 

for  p-  const,  1  -  >ty,  ^  and  with  determination  of  the  constant 

of  integration  and  constant  H  from  measurements  of  velocity  distri¬ 
bution  made  by  I.  Nlkuradse. 

To  obtain  a  universal  law  of  velocity  distribution  for  a 
turbulent  core  In  the  presence  of  gas  flow  with  heat  transfer.  It  Is 
necessary  when  Integrating  Eq.  (8)  to  take  Into  account  the  variation 
of  density  with  temperat\ire  and  to  Introduce  a  corresponding 
correction  for  heat  transfer  when  determining  the  constant  of  Integra¬ 
tion. 
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Integrating  £q.  (8),  taking  Into  account  variation  of  gas 

density  (l  »  xy,  t  ■  we  obtain 

I  .  I  r  du 

Our  measurements  Indicated  that  In  the  nea]>>wall  region  we  can 
assmne 


(10) 


Taking  Into  account  ( lO) ,  Eq.  (9)  reduces  to  the  form 


'  lny=_^  _L - - - iii  +  C  . 

**  l/^/.L  '  1-7’,  ^  i 


Taking  x  »  0.4  and  determining  the  constant  of  Integration  C 
by  Joining  the  turbulent  distribution  of  velocity  with  the  laminar 
distribution  of  velocity  in  the  Immediate  proximity  of  the  wall 
where  the  laminar  and  turbulent  tangential  stresses  are  of  the  same 
order  of  magnitude,  we  obtain _ 

j/  + 


ir-r. 


=  [5,7510^ Rey  j/^)  +  5,5] 


(11) 


,  ^  —Us.'  Re  — 

where  Tn'  “  •»  ' 


Constructing  velocity  profile  ( ll)  In  variables  ■  f(^)  for  a 
given  T,  we  obtain  a  network  of  curves  for  fixed  values  of  0^. 

If  the  experimental  velocity  values  In  the  test  cross  section  at  a 
corresponding  value  of  are  plotted  on  such  a  graph,  then, 
according  to  the  coincidence  of  the  velocity  profile  In  Its  logarithmic 
portion  with  one  of  the  curves  of  the  network,  we  will  arrive  at  a 
corresponding  value  of  c^.^  which  will  also  be  a  local  value  of  the 
friction  factor  at  the  wall  In  the  cross  section  under  consideration. 
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Having  completed  such  a  construction  for  the  various  form  para¬ 


meters  of  pressure  gradient 


Pi«i 


dx 


(12) 


If  8 

(where  —  v,-  is  the  kinematic  viscosity  of  the  potential 

flow  In  the  cross  section  under  consideration) ,  we  were  convinced 
that  the  logarithmic  law  of  velocity  distribution  Is  preserved  over 
approximately  one-third  of  the  boundary- layer  thickness  regardless 
of  the  pressure  gradient. 

An  analogous  construction  was  made  for  the  temperature  determina¬ 
tion.  A  universal  law  of  temperature  distribution  Is  preserved  over 
a  significantly  greater  portion  of  the  thickness  of  the  boundary 
layer  than  Is  the  vinlversal  law  of  velocity  distribution. 

The  second  method  of  determining  specific  heat  flow  was  the 
method  of  thermal  equlllbrliua. 

In  order  to  determine  the  distribution  of  local  values  of  tan¬ 
gential  stress  and  specific  heat  flow  over  the  cross  section  of  the 
boundary  layer,  the  differential  equations  of  momentum  and  energy  were 
Integrated,  taking  Into  account  the  equation  of  continuity  up  to  the 
current  value  of  the  y- coordinate. 

mie  momentum  equation  Is 


P« 


du 


du 


^  ,  I  dlrr) 
dx'r"  by  ’ 


the  equation  of  continuity  Is 


(15) 


d(rpu)  d(rpv)  _ 
dx 


and  the  energy  equation  Is 


97*,  ,  ar, 
p“^  +  P‘'aF 


SO.r 


d{rq) 

dy 


(1^) 

(15) 


where  x  and  £  are  coordinates  directed  along  ahd  normal  to  the  flow 
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respectively;  u  and  v  are  the  normal  and  tangential  components  of 
velocity.  After  Integrating  Eqs.  ( 13) ,  ( l4)  and  ( 15)  Integral  equa¬ 
tions  of  momentum  and  energy  were  obtained  for  a  boundary  layer  of 


thickness  from  0  to 

The  momentum  equation  Is 


+  du,  1  dr  I  dp,\ 

dx  \  w,  '  dx"^  r  "dx  Pi'  dx)^»~ 


ui-u  h(rpu)  T, 

dx  ‘^y 


The  energy  equation  Is 


Pi  Ml 


•(' 


ycosp\  T- 
■  '  'Pitt** 


(16) 


rf<p»  [1  du,  1 
dx  ■^17/57  +  7 


dr  1  rfp, 

Pi 'dx  '^Toi—T, 


dx 


1  d(Toi-TJi 


dx 


—  7*0  C  d  f  vj  1 

■  ^Pi  Mx  (foi  -  rj  J  dx  ‘^y  =■  8C,PiUi{T,i. 

0 


-TJ 


(17) 


loss  for  a 


where  0  =  1  I  (l  — ^^-^)rfythe  thickness  of  momentxim _ _ 

fixed  value  of  £;  and  the  displacement 

thickness  for  a  fixed  value  of 

Equation  ( 16)  Is  used  to  calculate  the  distribution  of  friction 
factor  c^  over  the  cross  section  of  the  boundary  layer  with  the  substi¬ 
tution  of  experimental  values  of  the  magnitudes  which  enter  Into  It. 
Heat  flows  for  a  fixed  value  of  ^  In  the  boundary  layer  are 
calculated  from  Eq.  ( 17)  . 

According  to  the  calculated  distributions  of  the  friction  factors 
and  heat  flows  In  the  cross  section  of  the  boundary  layer,  the 
corresponding  values  of  the  coefficients  of  turbulent  momentum  and 
heat  transfer  A^  and  A  were  determined  from  the  equations 


dy 


(18) 

(19) 
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5.  Results  of  the  Experimental  Investigations 

The  measurements  of  temperature,  total  and  static  pressure,  and 
heat  flow  In  five  Insulated  sections  of  the  experimental  diffuser  were 
made  under  22  regimes  of  motion  of  heated  air  and  In  4  regimes  with 
unheated  air.  Each  regime  Is  characterized  by  a  definite  air  tempera¬ 
ture  and  velocity  at  the  entrance  profile  of  the  test  section.  A 
change  of  entrance  air  parameters  by  regimes  made  It  possible  to 
significantly  alter  the  flow  rate  of  air  and  the  Reynolds  number. 

The  experiments  encompassed  Reynolds  niunbers  ( 1.688  +  8.48)  10®  and 
air  flows  of  0.35  to  1.37  kg /sec. 

The  presence  of  well-insulated  flow  stabilization  sections  and 
also  of  two  metal  nets  In  front  of  test  section  assured  equalization 
of  the  velocity  and  temperature  fields  over  the  flow  cross  section. 

The  origin  of  growth  of  the  boundary  layer  in  all  experimental  regimes 
was  related  to  the  entrance  cross  section  of  the  coOled  diffuser. 

6.  Analysis  of  Experimental  Data  ■ 

Figures  (3)  and  (4)  show  the  temperature  and  velocity  distribu¬ 
tion  over  the  cross  section  of  the  boundary  layer  for  various  values 
of  the  form  parameter  F.  It  is  apparent  fi?om  these  graphs  that  the 
velocity  profiles  In  the  diffuser  become  less  filled  with  an 
Increase  In  pressure  gradient.  The  same  tendency,  though  to  a  signi¬ 
ficantly  smaller  degree.  Is  observed  In  individual  temperature 
determinations.  Prom  a  comparison  of  the  graphs  of  velocity  and 
temperature  distribution  It  follows  that  In  gas  flow  with  a  positive 
p]*es8ure  gradient,  there  Is  no  similarity  of  velocities  and 
temperatures.  An  especially  significant  disruption  of  similarity 
of  these  fields  Is  observed  In  the  pre- separation  region. 
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Pig.  5.  Distribution  of  velocity  over 
the  cross  section  of  the  boundary  layer. 


Pig.  4.  Distribution  of  temperature  over 
the  cross  section  of  the  boundary  layer. 


6* 

Pigure  5  shows  the  dependence  of  the  form  parameter  H  ■  T" 

(with  calculation  of  9  in  the  variables  of  A.  Dorodnltsyn)  on  r. 

This  graph  shows  that  H  uniquely  depends  on  P,  where  upon  the  function 
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H  -  r(r)  is  Iniependent  or  the  Reynolds  number.  Therefore,  the  form 
parameter  r  expressed  by  E,.(l2)  l.  a  useful  form  parameter  for  the 
atudy  of  tbe  behavior  of  the  turbulent  boundary  layer  in  the  presenoe 
Of  a  flow  gradient  with  heat  transfer. 

The  connection  between  the  quantities  H  and  r  is  expressed  by 
equation 


Fig.  5.  H  vs.  f(-r) . 
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It  is  important  to  note  that  boundary- layer  separation  is  not 
observed  up  to  form  parameter  values  ‘greater  than  twice  the  separa¬ 
tion  values  of  r  ■  0.06  according  to  A.  Buri.  Accordingly,  the  values 
of  H  obtained  in  experiments  reach  2.5j  which  approaches  the  separa¬ 
tion  values  of  H  recommended  by  A.  E.  Doenhoff  and  N.  Tetervln  [8]  . 

The  actual  protraction  of  the  separation  point  in  comparison  with  the 
motion  of  uncompressed  fluids  is  explained  by  a  number  of  causes: 
the  stabilizing  Influence  of  the  cooled  walls  of  the  diffuser,  the 
axlsymmetrlcal  flow  under  our  conditions  and  also  by  the  fact  that 
in  the  experiments  of  I.  Nikuradse  [5],  the  results  of  which  were 
the  basis  of  Burl's  determination  of  the  separation  value  of  F,  the 
boundary  layer  fills  the  entire  cross  section  of  the  channel,  whereas 
our  experiments  were  conducted  with  a  potential  flow  over  a  signifi¬ 
cant  portion  of  the  passage. 

The  scattering  of  experimental  points  on  the  graph  of  H  »  f(r) 
may  be  explained  by  the  Influence  of  different  cooling  conditions. 
Apparently,  the  determination  of  H  with  respect  to  the  thickness  of 
momentum  loss  in  the  variables  of  A.  Dorodnltsyn  does  not  permit 
total  elimination  of  the  heat  transfer  effect. 

In  order  to  obtain  a  closed  system  of  equations  which  establish 
a  connection  between  the  characteristics  of  the  dynamic  boundary  layer, 
we  need  an  additional  form  parameter  which  defines  a  law  of  friction. 

On  the  basis  of  the  experimental  data  of  I.  Nikuradse  [5]*  A.  Burl 
showed  that  with  a  flow  gradient  of  an  uncompressed  fluid  without 
heat  transfer,  the  law  of  friction  may  be  expressed  by  the  equation 


;= 


(21) 


since  under  these  conditions  C  Is  a  single-valued  function  of  the  form 
parameter  F,  this  function,  like  the  function  H(F),  is  Independent  of 
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the  Reynolds  number. 

When  we  generalized  our  experimental  data  on  resistance  It 
became  clear  that  the  quantity  d  expressed  by  Eq.  (2l)  may  be  used  to 
establish  a  law  of  friction  also  In  the  presence  of  a  gas  flow  gradient 
with  heat  transfer  If  the  temperature  factor  Is  Introduced  Into 
the  functional  dependence  C  ■  f(r) . 

In  the  first  approximation  the  Influence  of  heat  transfer  on 
C  was  calculated  according  to  the  limit  formula  obtained  In  the 
literature  [ l6( sic) ] . 

Figure  6  shows  the  dependence  g  T  |  on  P  with  determination 

of  5  and  r  by  Eqs.  (2l)  and  (12)  .  The  significant  scattering  of 
experimental  points  on  this  graph  Is  apparently  explained  by  the  fact 
that  the  Indicated  limit  formula  was  obtained  for  a  nongradient  gas 
flow. 

The  experimental  points  in  Pig.  6  are  satisfactorily  distributed 
along  the  curve  of  the  equation 

5  =0,0128(1 

Taking  Eq.  ( 2l)  Into  account,  we  obtain  a  law  of  resistance 
T, _ C/._  0,0128  /  2  y  1 

Pi «?  “  2  ”  I  +  1/  (» -  0*®  •  (  23) 

Experimental  data  encompassing  a  variation  of  the  temperature 
factor  from  0.5  to  1,  and  values  of  form  parameter  F  from  0  to  0.12 
was  used  In  plotting  the  graph  In  Fig.  6. 

Figure  7  shows  the  dependence  of  form  parameter  H  on  the  form 
parameter  tj  Introduced  by  Qruschwltz  [2]: 

.  (24) 

which  Is  associated  with  the  form  parameter  ^7  the  equation 
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l*/« 


(25) 


obtained  according  to  the  exponential  law  of  velocity  distribution 

u  /  y  \ 

^r(v  • 

where  u(d)  Is  the  velocity  within  the  boundary  layer  at  a  distance 
y  ■  0  from  the  wall; 
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Pig.  6.  4  vs.  f(-r) . 


In  addition  to  our  own,  experimental  points  from  Oruschwltz  [2], 
Nlkuradse  [5],  and  Kehl  [6]  are  plotted  In  Pig,  7.  It  Is  apparent 
from  the  graph  that  the  experimental  points  of  Oruschwltz,  Nlkuradse 
and  Kehl  lie  nicely  on  the  Oruschwltz  curve  (Eq.  25)  In  the  region 
of  small  values  of  the  form  parameter  H,  Our  experimental  points 
cover  a  range  of  larger  values  of  the  form  parameters  H  and  t;.  In 
which  we  observe  a  significant  deviation  of  experimental  points  from 
the  Oruschwltz  dependence. 

Figure  8  shows  the  variation  along  the  longitudinal  coordinate 
of  local  values  of  the  resistance  coefficient  c^  at  the  wall  for  one 
of  the  regimes  (regime  13)  of  our  experiments.  The  curves  of  variation 
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of  the  resistance  coefficients  for  the  corresponding  conditions  based 
on  the  methods  of  Burl  and  Ye.  Kallkhman  are  also  plotted.  As  can 
be  seen  from  the  graph,  the  greatest  deviation  from  experimental 
points  Is  produced  by  the  Kallkhman  method  which  Is  based  on  the 
analogy  between  processes  of  heat  and  momentvun  transfer  In  a  boundary 
layer. 


Fig.  7.  H  vs.  f(T)) 
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7.  dynamic  Turbulent  Boundary  Layer  Calculation. 


Let  us  represent  Eq.  ( l)  in  the  following  form: 


Rei  'rf«i  ,  1  dr 


‘**''^*  I  tl4  \  1\  I  * 

-^  +  7 


‘r„Rei. 

3  • 

?\Ui 


where  Rcl=>^~^;L  ■  characteristic  length  of  measure;  and  m«o  " 

>  coefficient  of  dynamic  viscosity  of  the  potential  flow  at  stagnation 
temperature . 

Taking  into  account  Eq.  (2l)«  Eq.  (26)  is  reduced  to  the  form 

+  +  (27) 

We  will  designate  /'(r)=  — {^+ 1  pj. 

Figure  9  shows  the  dependence  of  the  function  P(r)  on  form  para¬ 
meter  r.  The  values  of  F(r)  calculated  from  the  values  of  C*  H, 
and  r  obtained  in  our  experiments  are  nicely  arranged  along  the 
line 

F  (0  =  0,0155- 3,55  r. 

ThuSj  although  the  quantities  C  and  H  substantially  depend  on  heat 
transfer  the  latter,  as  was  lobe  expected, iBs  practically  no  effect  on  the  func¬ 
tion  p(r) . 

Taking  into  account  Eq.  (29)  and  substituting  F  by  means  of  its 

expression  (12),  Eq.  (27)  is  transformed  into  a  differential  equation 

r/ii 

n-P  i-Vio  PI  r>a+-  rt-nHa-n  wt  l-Vi  no oi-ianp  RaiV  Integrating  thlS  eqUatlOn 


of  the  first  order  with  respect  to  Reg^' 
and  solving  for  Reg  we  obtain 


. .i-.F+c)]:” 


where  C,  the  constant  of  integration  is 


L  OiOjfM  Jj-;/ 


and  Xjj  is  the  value  of  the  dimensionless  coordinate  x  of  the  starting 
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point  of  the 


21 


On  the  basis  of  the  data  obtained,  Ihelr  treatment^  and  generaliza¬ 
tion,  which  made  It  possible  to  Integrate  the  equation  of  motion,  we  can 
recommend  the  following  sequence  of  calculation  of  a  dynamic  boundary 
layer  with  turbulent  gas  flow  In  an  axlsymmetrlcal  channel  with  a 
longitudinal  pressiire  gradient  and  with  heat  transfer: 

1.  With  given  laws  of  velocity  and  temperature  variation  of 
the  potential  flow  and  a  given  channel  radius  along  the  longitudinal 
coordinate  x,  the  Reg  number,  established  over  the  thickness  of 
momentum  loss  Is  determined  according  to  Eq.  (30)  as  a  function  of  x 
and,  consequently,  a  dependence  of  momentum  loss  thickness  d  as  a 
function  of  x  Is  determined. 

2.  The  magnitude  of  the  form  parameter  F  as  a  fxmctlon  of  x 
Is  determined  according  to  Eq.  (12). 

3.  The  local  value  of  the  resistance  coefficient  at  the  wall 

as  a  function  of  x  Is  determined  according  to  Eq.  (23)  For  a  given  law 
of  wall  temperature  with  respect  to  the  coordinate  x  and  for  a  known 
temperature  recovery  factor  for  the  potential  flow. 

4.  The  magnitude  of  the  form  parameter  H  as  a  function  of  x 
Is  determined  according  to  Eq,  (20). 

This  method  of  calculation  must  be  used  In  a  region  of  relatively 
small  temperature  factors  (?^  ■  0.5  to  l)  .  The  question  of  the 
Influence  of  more  Intense  cooling  or  heating  on  the  characteristics 
of  the  turbulent  boundary  layer  will  be  studied  In  the  future. 

8.  Thermal  Boundary  Layer  Calculation 

The  thermal  boundary  layer  Is  calculated  on  the  basis  of  an 
Integral  equation  of  energy  and  on  experimental  data  concerning  the 
distribution  of  temperature  In  the  boundary  layer  and  the  variation 
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of  the  temperature  of  the  channel  walls. 

Equation  ( 2)  Is  easily  reduced  to  the  form 


where 


£ 

llx 


x(rRe*''')  =  |-/'Re^St  Re^*, 


=  St  = 


.’I 


Pl'^l  {T 01  —  fa,) 


(31) 


Treatment  of  our  experimental  data  on  the  temperature  distribu¬ 
tion  Indicated  that  the  values  of  the  product  St  Re^'^^,  laid  out  as 
ordinates  for  the  values  of  the  form  parameter  of  the  thermal 
boundary  layer  with  a  pressure  gradient  equal  to 


(32) 


are  well  distributed  on  a  line  parallel  to  the  abscissa  (Fig.  lO) , 
thus  Indicating  an  Independence  of  the  product  St  Re^'^^  on  the 
pressure  gradient.  The  average  value  of  St  Re^'^^  In  the  graph  of 
Pig.  10  may  be  taken  as  equal  to  0.01. 
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Pig.  10.  Cij  vs.  f(r^  . 

Substituting  the  value  0.01  for  St  In  Eq.  (jl)  we  obtain 

1  ^  1  d  (fpi  — fg,) 

4  r  dx  4-  foi — T „  dx 

x(rRe‘'‘)-0,0126r  Re^,. 

Integrating  Eq.  (53)  with  respect  to  rRe^^^  and  then  expressing 

the  resulting  Integral  with  respect  to  Re^p  we  find 
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The  Stanton  number  as  a  functional  x  is  determined  according  to 


the  equation 


St  =  M. 
Roy* 


(55) 


The  thermal  boundary  layer  for  turbulent  gas  flow  In  an 
axl symmetrical  channel  with  a  longitudinal  pressure  gradient  Is 
calculated  In  the  following  sequence: 

1.  With  given  laws  of  velocity  and  temperature  variation  In  a 
potential  flow  as  well  as  a  given  wall  temperature  and  channel  radius 
along  the  longitudinal  coordinate  x,  Re^pCx)  and  then  the  dependence 
of  the  energy  loss  thickness  ^(x)  are  determined  by  Eq.  (34)  .  In 
the  case  when  no  analytical  law  Is  given  for  velocity  and  temperature 
variation  In  the  potential  stream  as  well  as  for  wall  temperature 
variation,  Eq.  (34)  may  be  Integrated  numerically. 

2.  Values  of  the  form  parameter  as  a  function  of  x  are  found 
from  Eq.  (32)  . 

3.  Values  for  Stanton  n\amber  are  calculated  by  Eq.(35)  . 

4.  Thermal  flow  at  the  wall  of  the  channel  Is  calculated  by 
equation 


~TJ St. 


5.  The  relation  between  the  form  parameters  F  and  F^  Is 
establJshed  on  the  basis  of  Eqs.  ( 12)  and  (32). 
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9.  Reynolds^ Analogy  on  Heat  and  Momentum  Transfer 


The  experimental  data  obtained  makes  it  possible  to  analyze  the 
question  of  the  accuracy  of  the  analogy  in  the  processes  of  heat  and 
momentum  transfer  in  the  presence  of  a  pressure  gradient. 

The  original  premises  for  Reynolds'  analogy,  as  is  well-known, 
are  the  equality  of  the  coefficients  of  turbulent  heat  and  momentum 
transfer  as  well  as  the  identity  of  the  laws  of  variation  of  t  and 
^  over  the  cross  section  of  the  boundary  layer.  The  result  of  these 
premises  is  the  similarity  of  the  velocity  and  temperature  profiles 
in  the  boundary  layer. 

A  significant  breach  in  the  similarity  of  the  velocity  and  temper¬ 
ature  profiles  in  the  boundary  layer  with  a  posit  Ive  pressure  gradient 
is  pointed  out  above.  In  order  to  clarify  the  causes  of  this  breach 
it  is  necessary  to  varlfy  the  accuracy  of  the  original  premises  of 
Reynolds'  analogy  under  conditions  of  a  gas  flow  gradient.  Figure  11 
shows  the  distribution  of  t  and  ^  over  the  cross  section  of  the  boundary 
layer  obtained  by  the  calculation  method  described  above.  As  can  be 
seen  from  the  graph,  the  laws  of  variation  of  t  and  over  the  cross 
section  of  the  boundary  layer  in  the  presence  of  a  gas  flow  gradient 
significantly  differ,  the  difference  being  greater  the  nearer  flow 
conditions  are  to  separation.  A  characteristic  maximum  with  a  sub¬ 
sequent  drop  to  zero  at  the  outer  edge  of  the  boundary  layer  is  ob¬ 
served  in  the  distribution  of  t  in  proportion  to  the  distance  from 
the  wall.  The  magnitude  of  specific  thermal  flow  drops  continuously 
from  maximum  at  the  wall  to  zero  at  the  outer  edge  of  the  boundary 
layer.  Thus,  one  of  the  premises  of  Reynolds'  analogy  is  not  verified 
by  the  experiments  with  a  gas  flow  gradient. 
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Pig.  11.  Distribution  of  heat  flow 
and  frictional  stress  over  the  cross 
section  of  the  boundary  layer  (regime 
12i  cross  section  2)  . 
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Pig.  12.  Distribution  of  the  cross 

section  boxmdary  layer  (regime  12,  cross 
section  2) . 

Pigure  12  shows  the  variation  of  the  relationship  of  the  coeffi¬ 
cients  of  turbulent  momentum  and  heat  transfer  ( turbulent  Prandtl 
number)  over  the  cross  section  of  the  boundary  layer.  Prom  Pig.  12 
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it  follows  that  near  the  wall,  the  coefficients  and  A,  are  the 

q  T 

same.  The  coefficient  of  turbulent  heat  transfer  exceeds  the  coeffi¬ 
cient  of  momentum  transfer  In  proportion  to  the  distance  from  the  wall, 
and  at  some  distance  from  the  wall  A^  Is  approximately  twice  A^. 
Subsequently  the  ratio  A^/A^  drops  In  proportion  to  distance  from 
the  wall.  Consequently  the  second  premise  of  Reynolds' analogy  (on 
the  equality  of  the  coefficients  A^  and  A^)  Is  also  not  donflrmed 
by  experimental  data. 

The  regularity  of  the  variation  of  A^/A^  over  the  cross  section 
of  the  boundary  layer  which  we  obtained.  Is  In  qualitative  agreement 
with  experimental  data  of  Ludwlng  [9],  Elias  [10],  and  Peydzh  and  Pokner 
[11].  In  particular,  the  latter  obtained  A^A^  =  2  on  the  basis  of 
the  measurement  of  temperature  and  velocity  distribution  In  a  free 
flow.  It  may  be  assumed  that  at  some  distance  from  the  wall  the 
mechanism  of  heat  and  momentm  transfer  obeys  the  regularities  of 
free  turbulence.  Prom  the  physical  point  of  view  the  Influence  of 
free  turbulence  must  predominate  In  proportion  to  the  nearness  of  the 
points  of  separation,  since  in  this  case  the  influence  of  the  wall 
on  the  characteristic  of  turbulence  In  the  boundary  layer  decreases. 

The  analysis  which  was  carried  out  Indicates  that  the  Reynolds 
analogy  is  not  verified  for  gas  motion  with  a  longitudinal  pressure 
gradient.  Consequently  all  of  the  analytical  methods  of  calculating 
the  turbulent  boundary  layer  which  are  based  on  the.  use  of  Reynolds' 
analogy,  Kallkhman's  In  particular,  cannot  be  extended  to  the  case  of 
gas  motion  with  a  significant  positive  pressure  gradient. 

Conclusions 

1.  As  a  result  of  the  experimental  studies  carried  out,  new 
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experimental  data  were  obtained  on  the  characteristics  of  a  tvirbulent 
bovindary  layer  with  gas  motion  In  an  auclsymmetrlcal  diffuser  with 
cooled  walls. 

2.  Analysis  of  the  experimental  data  obtained  Indicated  that 
existing  methods  of  turbulent  boundary  layer  calculation  (of  Burl^ 
Gruschwltz  and  Kallkhman)  are  not  verified  by  experiments  with  gas 
motion  In  an  axlsymmetrlcal  diffuser  with  cooled  walls. 

3.  Experimental  values  of  the  form  parameter  r  corresponding  to 
the  separation  point  of  a  turbulent  boundary  layer  significantly 
exceed  the  corresponding  values  according  to  the  data  of  I. 

Nlkuradse . 

4.  Experimental  values  of  the  form  parameter  C  He  between 

the  corresponding  values  obtained  by  the  methods  of  Burl  and  Gruschwltz. 

4 .  Experimental  data  were  obtained  on  the  dependence  of  the  form 
parameter  H  on  the  Gruschwltz  form  parameter  n  in  a  region  which  was 
not  previously  covered  by  experimental  data. 

6.  Analysis  of  experimental  data  Indicated  that  the  Burl  method 
may  be  extended  to  the  case  of  gas  motion  in  an  axlsymmetrlcal  diffuser 
with  cooled  walls  by  the  introduction  of  suitable  refinements  for 

the  graphs  of  H  *  f(r) ;  ^  »  f(r)and  P(r)»  f(r)  . 

7.  The  proposed  method  of  calculation  is  found  to  be  In  good 
agreement  with  experimental  data  both  with  regard  to  the  calculation 
of  the  Integral  characteristics  of  the  turbulent  boundary  layer  auid 
to  local  values  of  the  coefficient  of  friction. 

8.  Of  the  existing  methods  of  turbulent  boundary  layer  calcu¬ 
lation,  the  method  based  on  the  premises  of  Reynolds  analogy  (the 
method  of  L.  Ye.  Kallkhman)  gives  the  greatest  divergence  from  experi¬ 
ments. 
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9.  Experiments  Indicated  that  with  a  gradient  of  gas  motion  there 
Is  observed  a  significant  breach  of  similarity  between  the  velocity 
and  temperature  profiles. 

10.  The  ratio  q/r  cannot  be  taken  as  constant  over  the  cross 
section  of  the  bOTindary  layer  with  a  gradient  'of  gas  motion. 

11.  The  ratio  of  the  coefficients  of  turbulent  heat  and  momentum 

transfer  varies  significantly  over  the  cross  section  of  the 

boundary  layer  from  A^/A^  ■  1  at  the  wall  to  A^/A^  ■  2  In  the  main 
section  of  the  turbulent  boundary  layer. 
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